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ABSTRACT 

Using a resimulation technique, we perform high-resolution cosmological simulations of dry mergers 
in a massive (1O 15 M ) galaxy cluster identified in the Millennium Run. Our initial conditions include 
well resolved compound galaxy models consisting of dark matter halos and stellar bulges that are 
used to replace the most massive cluster progenitor halos at redshift z = 3, allowing us to follow the 
subsequent dry merger processes that build up the cluster galaxies in a self-consistent cosmological 
setting. By construction, our galaxy models obey the stellar mass-size relation initially. Also, we study 
both galaxy models with adiabatically contracted and uncompressed halos. We demonstrate that the 
brightest cluster galaxy (BCG) evolves away from the Kormendy relation as defined by the smaller 
mass galaxies (i.e., the relation bends). This is accompanied by a significantly faster dark matter mass 
growth within the half-light radius of the BCG compared to the increase in the stellar mass inside 
the same radius. As a result of the comparatively large number of mergers the BCG experiences, its 
total mass-to-light ratio becomes significantly higher than in typical elliptical galaxies. We also show 
that the mixing processes between dark matter and stars lead to a small but numerically robust tilt 
in the fundamental plane and that the BCG lies on the tilted plane. Our model is consistent with 
the observed steepening of the logarithmic mass-to-light gradient as a function of the stellar mass. 
As we have not included effects from gas dynamics or star formation, these trends are exclusively 
due to TV-body and stellar-dynamical effects. Surprisingly, we find only tentative weak distortion in 
the Faber- Jackson relation that depends on the aperture size, unlike expected based on studies of 
isolated merger simulations. This may be due to differences in the distribution of galaxy orbits, which 
is given in our approach directly by the cosmological context while it has to be assumed in isolated 
merger simulations, and the fact that the BCG is located deep in the cluster potential well. Another 
uncertainty in both approaches lies in the definition of the spatial extent of the BCG. 

Subject headings: clusters: general - galaxies: elliptical 



1. INTRODUCTION 
1.1. What makes BCGs special ? 

The brightest cluster galaxies (BCGs) are special. 
They are the most massive and luminous galaxies in the 
Universe. They are typically located in the very centers 
of clusters of galaxies which indicates that their forma- 
tion is closely linked to that of the clusters themselves. 
Their formation history is therefore distinct from typical 
elliptical galaxies (e.g., Lin & Mohr 2004, Brough et al. 
2005, De Lucia & Blaizot 2006 and references therein). 
Recent analysis by Best et al. (2007) shows that BCGs 
are also more likely to host active galactic nuclei than 
other galaxies of the same stellar mass. This indicates 
that these objects play a pivotal role in quenching cooling 
flows and star formation in clusters. 

BCGs have been recently shown to lie off the standard 
scaling relations of early-type galaxies (Lauer et al. 2007, 
von der Linden et al. 2007, Bernardi et al. 2007). In 
particular, they show excess luminosity (or stellar mass) 
above the prediction of the standard Faber- Jackson rela- 
tion at high galaxy masses. This hints at the interesting 
possibility that their black hole masses may be larger 
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1.2. Formation of BCGs in the cosmological context 

Early theoretical studies suggested that BCGs are 
formed as a result of star formation in cooling flows 
(Fabian 1994). However, recent X-ray observations per- 
formed with Chandra and XMM Newton show that cool- 
ing rates in nearby clusters are too low to explain the 
masses of BCGs. Tidal stripping and dynamical fric- 
tion acting on the cluster galaxies have been suggested as 
an alternative BCG formation mechanism (Oemler 1976, 
Richstone 1976) but this mechanism turns out to be too 
slow due to the low cross section for galaxy-galaxy inter- 
actions in virialized cluster centers where relative veloc- 
ities of galaxies are large (Ostriker 1980). Nevertheless, 
this cannibalism may be responsible for the formation 
of diffuse stellar envelopes extending up to hundreds of 
kiloparsecs (the 'cD' galaxies). 

However, in a hierarchical structure formation theory, 
groups of galaxies will form before the formation and viri- 
alization of a cluster. Such groups have sufficiently small 
velocity dispersions that a significant number of galaxy- 
galaxy mergers occurs. This process should preferentially 
lead to the formation of elliptical galaxies in high galaxy 
density environments. Stellar population synthesis mod- 
els demonstrate that the bulk of star formation in most 
massive galaxies took place prior to z ~ 2 (Thomas et 
al. 2005, Treu et al. 2005, Jimenez et al. 2006). Semi- 
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analytical models also show that stars that make up most 
of the BCG mass are formed very early on (80% by z ~ 3; 
De Lucia et al. 2006, De Lucia & Blaizot 2007). Only af- 
ter most stars have been formed docs the final galaxy as- 
sembly take place (the downsizing phenomenon; Bundy 
et al. 2005, Pannella et al. 2006, Calura et al. 2008, 
but see the revised monolithic collapse model of Chiosi 
& Carraro (2002) and the results of Cimatti et al. 2006 
who, based on the same COMBO-17 and DEEP2 data as 
Bell et al. 2004, argue for "top-down" assembly of early- 
type galaxies, i.e., that downsizing concept may have to 
be extended to the mass assembly itself as the build-up 
of the most massive galaxies precedes that of the less 
massive ones.). Therefore, the final galaxy mergers in 
cluster formation are expected to be predominantly dis- 
sipationless (or dry, i.e. not involving large amounts of 
gas; Khochfar & Burkcrt 2003). This picture is also con- 
sistent with the analysis based on COMBO-17 results of 
Bell et al. (2004) who show that the stellar mass on the 
red sequence increased passively by a factor of ~ 2 since 
z ~ 1. We note that recent theoretical results suggest 
that the amount of gas associated with major mergers in 
the early Universe was much larger then at the present 
epoch (Khochfar & Silk 2006). Such wet mergers led to 
massive starbursts early on. However, at lower redshifts 
the amount of gas was much lower and the dominant 
mechanism for the galaxy growth in mass and size were 
dry mergers (e.g., Hopkins et al. 2008). This is evidenced 
indirectly by the strong size evolution of the most mas- 
sive galaxies since z = 2 (Trujillo et al. 2007) while 
preserving the mass-size relation found by Trujillo et al. 
(2004). 

Despite these findings, we note that the detailed pic- 
ture of the giant elliptical assembly is further constrained 
by photometric and chemical observables. They suggest 
that it is difficult to explain giant elliptical by a pure se- 
quence of multiple minor dry mergers or via major dry 
mergers (i.e., some admixture of wet mergers may be re- 
quired to properly model chemical observables; Pipino 
& Matteucci 2008). However, the scatter in <r-[Mg/Fe] 
relation is large enough to permit one to three major dry 
mergers during the galactic lifetime (Bell et al. 2006). 
While the majority of BCGs from the sample studied by 
Bildfcll et al. (2008) have bluer colors with increasing 
radius, about 25% possess bluer cores. The existence of 
such cores indicates that some cold gas must have been 
replenished in that subsamplc of BCGs at late epochs via 
minor wet mergers. Interestingly, Bildfell et al. (2008) 
find that the Kormendy relation of the BCGs is steeper 
than that of the local ellipticals. Kaviraj et al. (2008) 
argue that bright early-types potentially form 10-15 per 
cent of their stellar mass since z = 1 and, according to 
Pipino et al. (2008), the recent star formation in the blue 
core BCGs typically has an age less than 0.5 Gyrs and 
contributes mass fractions of less than a percent, von 
der Linden et al. (2007) show that BCGs exhibit higher 
[a/Fe] values than same mass ellipticals and claim that 
this indicated that star formation may have occurred over 
a shorter timcscale in BCGs. This then suggests a possi- 
bility that the [a/Fe] radial gradients in BCGs could be 
consistent with being comparable to those in same mass 
ellipticals (Brough et al. 2005) rather then being unable 
to survive repeated dry mergers. 

We also point out that the AGN feedback plays an im- 



portant role in explaining the early and rapid formation 
of most massive ellipticals by shutting off star formation. 
The effect of AGNs may be felt more strongly in more 
massive galaxies that have been heated by virial shocks. 
This made their gas more dilute and more vulnerable 
to AGN feedback leading to passive evolution and 'red- 
and-dead' massive spheroids (Dekel & Birnboim 2006). 
Recent high-resolution FLASH code simulations of the 
interaction of rclativistic jets with the ISM confirm that 
the star formation rate will be inhibited although it can 
increase on short time-scales of the order of 10 5 — 10 6 yr 
Antonuccio-Delogu & Silk 2008). 

As mentioned above, recent analysis of the SDSS and 
Hubble data demonstrates that BCGs form a separate 
class of objects that depart from the usual Faber- Jackson 
relation. This effect has been seen in numerical simu- 
lations of dry mergers performed by Boylan-Kolchin et 
al. (2006) who studied collisions of isolated galaxy pairs. 
They demonstrated that the initial orbital parameters 
play an important role in shaping the Faber-Jackson re- 
lation. In their model, galaxies that collide on more ra- 
dial orbits, such as those present during the cluster as- 
sembly, lead to more puffed up merger remnants. Such 
final merger products are characterized by lower inter- 
nal velocity dispersions than those resulting from lower 
cllipticity collisions. This mechanism is consistent with 
the fact that most massive galaxies experienced at least 
one merger since z < 1.5 (Conselice et al. 2007). As 
noted by Lauer et al. (2007), the increase in the ratio of 
the velocity dispersion of galaxies moving in the cluster 
potential to the internal velocity dispersion of individual 
cluster galaxies may also contribute to the puffing up ef- 
fect. This is because some of the orbital energy of the col- 
liding galaxies will be transferred to the internal energy 
of the merger remnant. Another effect that may influ- 
ence the slope of the Faber-Jackson effect is tidal heating 
of galaxies moving in the cluster potential. It is due to 
all of the above complexities that one has to resort to 
cosmological simulations involving both the dark matter 
and stellar components to capture all relevant processes 
leading to the formation of BCGs. It is the purpose of 
this paper to address some of the above issues. We will 
deliberately use purely collisionlcss simulations, allowing 
us to obtain reliable A-body results in a fully cosmologi- 
cal context without having to delve into the complexities 
of gas dynamics and star formation. However, we note 
that the approach we outline here can be used also in fu- 
ture hydrodynamic simulations of the cluster formation 
process. 

The paper is organized as follows. In Section 2, we out- 
line the numerical approach used to simulate the cluster 
and its galaxies. Section 3 discusses our results (scaling 
relations, tilt in the fundamental plane, the evolution 
of the mass-to- light ratio, and merger rates). Then in 
Section 4, we discuss the numerical convergence and the 
dependence of our results on the initial conditions and 
the galaxy models. Finally, Section 5 summarizes our 
findings. 

2. METHODS 

Our approach to simulate the evolution of the BCG 
is similar in spirit to that of Dubinski (1998), who re- 
placed 100 dark halos at z = 2 in a proto-cluster region 
with galaxy models constructed in isolation, and then 
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continued the simulation to the present epoch. However 
our method differs in several important details, and for 
dcfinitcness, we list its primary steps in the following: 

• We identify a massive 10 15 M Q cluster at z = in 
the Millennium Run (Springel et al. 2005a) and 
track its particles back to the initial conditions. 

• We determine the Lagrangian region of the cluster 
material in the original unperturbed initial condi- 
tions and produce a new refined unperturbed par- 
ticle load at z = 127 which is adequate for a zoom 
simulation onto the selected target halo at higher 
resolution. 

• We generate the actual initial conditions by per- 
turbing the particle load obtained in the previous 
point with the Millennium displacement field, aug- 
mented with additional small-scale waves that can 
now be added due to the improved mass resolution 
of ~ 10 9 h- l MQ/F^ ola . We consider F zoom = 2 
and F zoom = 4 in this paper, corresponding to ap- 
proximately 8 and 64 times the original resolution. 
Further away from the cluster, the resolution is 
progressively downgraded in mass, and the grav- 
itational softening lengths are enlarged, ensuring 
however that the tidal forces acting on the high- 
resolution region are represented accurately at all 
times. 

• We evolve these initial conditions to z = 3.0. 

• We identify the fifty most massive dark matter ha- 
los at this redshift (almost all of them end up in the 
final cluster) and replace them with galaxy models 
at either the same mass resolution as in the high- 
resolution region, or with the resolution enhanced 
by another factor of eight. The models consist of 
both dark matter and stellar particles (see below). 
In every run the dark matter and star particles of 
galaxy models have the same mass in order to min- 
imize differential two-body heating effects between 
the two components. The maximum mass and force 
resolutions for each of our simulations are listed in 
Section 3. 

• We evolve the cluster from z = 3.0 to the present 
epoch. 

We note that this approach is also similar in spirit 
to that recently employed by Rudick, Mihos & McBride 
(2006) who focused on the problem of the intracluster 
light and studied significantly less massive (~ 10 14 M Q ) 
clusters. We only consider spheroidal galaxy models in 
this paper that include a stellar bulge and a dark matter 
halo. The bulge size obeys the size-stellar mass relation 
of Shen et al. (2003): 

/ M \ °- 56 

where R e is the half-light radius, and M sta r is the stellar 
mass of a given galaxy. There is observational evidence 
(Trujillo et al. 2004, Mcintosh et al. 2005) that this 
relation remains unchanged up to z ~ 1 as long as the 
most massive ellipticals are excluded. Note that more 



recent observations suggest that the size-luminosity rela- 
tion evolves with redshift for high-mass ellipticals (Zirm 
et al. 2007, van Dokkum et al. 2008) but at these mass 
scales this trend is naturally explained by our simulations 
as explained in section 3). For the sake of simplicity we 
here assume that the relation also holds at our redshift 
(z = 3) for constructing the initial galaxy models (see 
Section 3 for more discussion of this assumption). We 
describe the distribution of mass in stars m*(r) in the 
bulge by a Hernquist profile, i.e., 



where a is a characteristic radius, M» is the total mass 
of the stellar component in a given galaxy and r is the 
distance from its center. We truncate the Hernquist pro- 
file at r = 10 a in our numerical realizations in order 
to avoid the occurrence of a very small number of very 
distant star particles. For this choice of parameters, the 
half-light radius is R e = 1.339 a. 

We assume that 90% of the original dark matter mass 
of each halo excised from the zoomed Millennium sim- 
ulation consists of dark matter and the rest contributes 
to the stellar component. The dark matter halos are 
modeled with the NFW profile (Navarro, Frenk & White 
1996) truncated at r2oo, with a concentration parame- 
ter c v ir given by a modified version of the Bullock et al. 
(2001) prescription, 

c vir (/i, a) = K—, (3) 
a c 

where K is a constant, a = 1/(1 + z), a c is the epoch 
of collapse of a given halo and \x = M v - 1T (a)/M*(a) is a 
ratio of a halo mass and the typical halo mass at the 
same epoch a. The typical mass at an epoch a is de- 
fined as a{M*(a)] — 1.686/D(a), where D(a) is the lin- 
ear growth factor and a(M) is the rms density fluctua- 
tion on a comoving scale containing the mass M. The 
collapse epoch is defined by M*(a c ) = FM v ; r , where F 
is a free parameter. In the computation of c(/i,a) V ir, as 
well as in all of our simulation runs (carried out with the 
GADGET-3 code; Springel 2005c), we assume the cos- 
mological parameters of the Millennium Run, given by 
Q m = 0.25, Q\ = 0.75, a scale invariant slope of the 
power spectrum of primordial fluctuations (n = 1.0), a 
fluctuation normalization erg = 0.9, and a Hubble con- 
stant of H = lOO/ikms^Mpc- 1 = 73kms~ 1 Mpc -1 . 
The constants K — 3.4 and F = 0.01 have been cho- 
sen following Maccio et al. (2006). We note that re- 
cently Gao et al. (2008) have provided a yet more ac- 
curate parametrization of the concentration dependence 
on mass and redshift than given by equation ([3]), but our 
results are insensitive to the size of this revision. 

Once the density structure of our galaxy models 
has been specified, velocity distributions of the bulge 
and dark matter particles are obtained by solving the 
Jeans equations following the method of Springel et al. 
(2005b). For the sake of simplicity, we assume that the 
dark matter halos have zero net angular momentum. 

Where indicated, the dark matter distribution has 
been modified by including adiabatic contraction follow- 
ing the classic Blumenthal et al. (1986) approximation. 
Hydrodynamical simulations by Gnedin et al. (2004) sug- 
gest that this approximation overpredicts the amount of 
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contraction at small radii. We note however that these 
differences between hydrodynamical simulations and the 
contraction model are smaller than the differences be- 
tween the uncompressed and contracted halos, and that 
an overestimate of the contraction effect allows us to 
bracket reality between the two models we consider. 

Nevertheless, we want to stress that our toy galaxy 
models are too simplistic to represent real galaxies very 
accurately. However, the advantage of our approach is 
that it gives us better control over the parameters of the 
galaxy population under study and allows us to cleanly 
separate the pure A^-body effects from the complicated 
physics of gas dynamics, radiative cooling, star forma- 
tion, etc. Our method should still be sufficiently accurate 
to investigate the evolution of the main trends in galaxy 
scaling relation, especially in the if-band. While the lu- 
minosity of the early type galaxies on the red sequence 
does evolve with redshift, the passive stellar evolution 
does not change the if -band luminosity much and hence 
should not dramatically affect our results. This approach 
differs from, e.g., that of Conroy et al. (2006) and Berrier 
et al. (2008) who combine dark matter simulations with 
recipes for assigning light to dark matter halos and study 
the assembly of a large number of clusters, or from that of 
Murante et al. (2007) who incorporate the gasdynamical 
effect and star formation at the expense of lower dynam- 
ical range. 

In our analysis, the BCG was identified using first the 
friends-of-friends algorithm to find the cluster halo and 
its galaxies, and then by applying the SUBFIND group 
finding algorithms to look for gravitationally bound 
structures. We did not set any strict isophotal limit when 
defining galaxies, instead we let SUBFIND identify the 
bound stellar and dark matter part. This sets a certain 
isophotal limit for the BCG. The corresponding isopho- 
tal limit for non-BCGs may be slightly different but that 
should not affect their half-light radii significantly as they 
do not possess extended envelopes. Moreover, unless the 
change in the isophotal limits is huge, it is not likely to 
substantially affect the total luminosities and half-light 
radii of the galaxies in our sample. 

3. RESULTS 

All results presented in this section correspond to our 
default case where the initial galaxy models were adi- 
abatically contracted. For details of the numerical pa- 
rameters for all runs see Section 4. Figure 1 shows the 
evolution of the relation between the half-light radius R e 
and the stellar mass within that radius. This relation is 
similar to the Kormendy relation. The initial (z — 3) and 
final (z = 0) models are represented by the red and black 
points, respectively. Intermediate redshifts are marked 
by open circles with colors as given in the inset. The fig- 
ure demonstrates that, while smaller mass galaxies show 
relatively little evolution, the BCG clearly gains a sub- 
stantial amount of mass and increases in size. Interest- 
ingly, it clearly evolves off the extrapolation of the initial 
radius-mass relation. We stress that, in the current sce- 
nario, this effect is entirely due to A^-body process and 
is not caused by gas-dynamical processes, radiative cool- 
ing, star formation, etc. This effect has been reported 
in the observational work of Lauer et al. (2007). They 
showed that the half-light radii of BCGs lie above the 
log-log relation between the half-light radius and stellar 



luminosity defined by smaller mass galaxies. They also 
attribute this effect to a " progressive change in the char- 
acter of 'dry mergers' at higher galaxy masses". This 
has also been reported even more recently by Bildfell et 
al. 2008 who state that "the Kormendy relation of our 
BCGs is steeper than that of the local ellipticals, suggest- 
ing differences in the assembly history of these types of 
systems." We also note that the trend for the size of the 
most massive galaxies to evolve with redshift (Zirm et al. 
2007, van Dokkum et al. 2008) is consistent with the evo- 
lution of our simulated relation. Interestingly, our results 
also imply the steepening of the logarithmic mass-to- light 
gradient as a function of the stellar mass. The logarith- 
mic M/L gradient (Napolitano et al. 2005) evaluated 
using r- m = R c and r out = R V1I is approximately V;7 « 
(R e /R viT )[(M UM /M s ) out - (M DM /M s ) in ] cx i? e /i?vir be- 
cause R c < i? vir , (M DM /M s ) out > (M DM /M s ) in and 
(Mr>M/M s ) out is a constant (Mdm and M s denote the 
dark matter and stellar masses enclosed within a given 
radius, respectively). Since R c cx Mf 1 we have V/7 cx 
Mf 1//3 . Now, for a = 0.56 the relation is flat while 
for higher stellar masses a increases and so does the log- 
arithmic M/L gradient. This is in agreement with the 
observational findings of Napolitano et al. (2005). 

Figure 2 shows the evolution of the stellar (violet) and 
dark matter (red) components within the half-light ra- 
dius of the BCG as a function of redshift. The values 
on the vertical axis are normalized to the dark matter 
mass within the half-light radius at z = 3. It is evident 
from this figure that the bend in the R e — M relation is 
accompanied by a significant gain in mass of the BCG 
both in terms of the stellar and dark matter components. 
Interestingly, the ratio of the dark matter-to-stellar mass 
increases with time significantly. In other words, the re- 
gions of the BCG within the half-light radius become dis- 
proportionally more massive with time for a given stellar 
luminosity. We note that the stellar mass accretion rates 
of BCGs (a factor of ~ 2.5 since z ~ 2) obtained by 
Romeo et al. (2008) are consistent with our findings. 

The change in the dark matter-to-stellar mass ratio 
is likely caused by the large number of mergers that the 
BCG experiences over time, especially when compared to 
other galaxies. In Figure 3 we show a histogram of the 
number of mergers experienced by the BCG, the second 
largest and the third largest galaxy in the cluster. A 
merger is here defined as an event where over 50% of the 
stellar mass from one galaxy ends up in another galaxy. 
The vertical axis is for the cumulative number of mergers 
from the initial redshift z- m { = 3 down to a given redshift 
z < z in i (shown on the horizontal axis). As is clearly 
seen, the number of mergers declines rapidly with galaxy 
mass, and most of the mergers happening in the cluster 
are between the BCG and smaller galaxies. 

We note that the results for the most massive galaxies 
are quite insensitive to the exact definition of the n-th 
most massive galaxy, which can be defined as either (i) 
n-th largest at a given redshift or (ii) n-th largest at 
z = 0. Strictly speaking, we adopt the latter definition 
and backtrace the history of the particles belonging to 
such-defined n-th largest galaxies in the computational 
volume. The fact that the results for the most mas- 
sive objects do not depend on the choice of the defini- 
tion shows that the most massive galaxies preserve their 
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identities and it is mostly smaller objects that accrete 
onto them. The Figure 3 indeed confirms that the BCG 
experiences overwhelmingly more mergers than smaller 
galaxies. This fact is the key reason for the change in the 
mass ratio between dark matter and the stellar compo- 
nent within the BCG's half-light radius. The number of 
mergers experienced by the most massive galaxy is zero 
between z = 2 and z = 3 and the smaller the galaxy mass 
the lower the redshift at which it experiences a merger 
(if any). Therefore, the results are not sensitive to the 
exact insertion redshift for the galaxy models. 

We note that the boost in the dark matter to stellar 
mass ratio could be contaminated by mergers of small 
starless halos (which were not replaced with galaxy mod- 
els) with a bigger halo. However, the dynamical friction 
time is a strongly decreasing function of the mass of the 
object that is merging with a bigger one. Therefore, the 
contribution to the growth in the BCG mass that comes 
from the small halos that we do not populate with stars 
is expected to be weak. This fact is not inconsistent with 
the above observation that the most massive objects do 
not merge with the BCG. It simply means that among 
the objects within a certain mean distance from each 
other, the most massive ones are the first ones to merge. 
In other words, the biggest galaxies in our simulations 
have simply been prevented from merging thus far by 
their large distance. 

Curiously, the stellar mass-c dispersion shown in Fig- 
ure 4 does not show strong bend similar to that presented 
in Figure 1. Figure 4 shows the evolution of the stellar 
mass within the half-light radius and the velocity disper- 
sion a which, in our model, is essentially an equivalent 
of the Faber- Jackson relation. The initial slope of this 
relation is ~ 3.62. The meaning of the symbols is the 
same as in Figure 1. The stellar velocity dispersions are 
measured within half-light radii, are three-dimensional 
(i.e., all three components of velocity are used to com- 
pute them) and are intensity- weighted (i.e., in our case, 
weighted by the stellar mass integrated along the line-of- 
sight). At the final redshift the BCG lies on the relation 
defined by the extrapolation of the initial galaxy models 
to more massive objects. 

The right panel shows the flux- weighted a as a function 
of stellar mass for three different apertures: half-light 
radius (red points), 1/8 or the half-light radius (green 
points) and 1/16 of that radius (blue points). There is a 
systematic trend for the BCG to move to smaller veloc- 
ity dispersions as the size of the aperture is decreased. 
At the same time, low-cr objects only increase their scat- 
ter around the unaltered mean values. There is also a 
tentative weak hint of a bend in this relation but this 
conclusion depends on the precise measurement of low- 
est a galaxies that show scatter and the magnitude of 
the bend may be slightly affected by the precise defini- 
tion of the BCG boundary that separates it from the 
intracluster light (ICL). 

We also note that the dependence of sigma on luminos- 
ity is slightly more subtle. In the case of size-luminosity 
relation, for a given binding energy of the merging pair, 
the satellites on high angular orbits need to lose more an- 
gular momentum to merge and then end up being more 
compact. Those on lower angular momentum orbits are 
then more puffed up. On the other hand, the velocity dis- 
persion close to the galactic center is strongly influenced 



by the stellar bulge and declines with radius up to a point 
where the cluster potential starts to dominate leading to 
the increase of the dispersion with the distance from the 
BCG center. This means that the cluster potential may 
try to slow down the increase of the slope with luminos- 
ity of the Faber-Jackson relation. In any case, fitting the 
extended wings of the BCG light distribution, beyond of 
the extent of this central galaxy as identified in our on- 
the-fly algorithm, could perhaps increase the bend in this 
relation (and the Kormendy relation) slightly and allow 
one to study the evolution of the ICL but such effects 
should not dramatically alter our conclusions. We defer 
such detailed analysis to a forthcoming publication. 

Liu et al. (2008) analyzed a sample of BCGs and ellip- 
ticals and found that the slope of both the size- luminosity 
and the Faber-Jackson relation changes with the galaxy 
luminosity. However, they also analyzed BCGs and 
bright ellipticalls falling into the same magnitude range. 
Curiously, they found that, while "the power-law indices 
for the size-luminosity relations for BCGs are steeper 
then those for bright non-BCG elliptical galaxies, the 
power-law indices for the Faber-Jackson relations are also 
steeper, but only at ~ la levels, and so are not statis- 
tically significant" (sec their Figure 14). This trend ap- 
pears to be consistent with our simulations. The depar- 
ture of the BCG from the original log-log size-luminosity 
relation increases with time and the slope of the size- 
luminosity for the BCG also seems to rise. This clearly 
suggests that an ensemble of simulated clusters would 
show a different size-luminosity slope for BCG and non- 
BCG galaxies in the same luminosity range, but a reli- 
able determination of the size of the effect would require 
a much larger simulation set. 

Finally, Figure 5 shows the "fundamental plane" corre- 
sponding to the initial models at z = 3 (red points) and 
the galaxies at the final redshift z = (violet points). 
The "brightness" parameter is I e = m stms / (ttR^) : where 
w s tars i s the stellar mass within the half-light radius. The 
absolute values on the axes are not important here as we 
are only interested in the evolution of the slope of this 
relation. The initial slope of the relation is ~ 0.83 (by 
definition, the virial relation corresponds to 1.0). Indeed, 
there appears to be a small but numerically robust tilt 
in the fundamental plane. The final slope of the rela- 
tion is ~ 0.72. The BCG seems to lie on the tilted plane 
but, strictly speaking, this conclusion depends on the 
robustness of the values measured for the least massive 
galaxies (see below for a discussion of numerical con- 
vergence). Nevertheless, the trend for the most massive 
object(s) to move away from the original relation is real. 
We stress again that this effect is entirely due to stellar 
dynamical/iV-body effects. 

4. NUMERICAL CONVERGENCE AND MODEL 
ROBUSTNESS 

We performed several test runs in order to assess the 
numerical convergence of our results and their sensitivity 
to the initial galaxy models. For testing the basic numer- 
ical convergence we in general used uncontracted galaxy 
models; these have less particles inside the half-light radii 
and therefore yield more conservative convergence tests. 

Figure 6 shows one of the results of these tests. It 
presents the relation between the half-light radius and 
the stellar mass within that radius. Initial models arc 
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denoted by open red circles. All other points are for the 
galaxies at the final redshift z = 0. Different sets of 
points grouped by color correspond to different resolu- 
tions. In this figure, we show the results for four different 
runs: 

• run A: maximum mass resolution at the base- 
line level of the zoom simulation equal to 1.25 x 
10 8 h~ 1 M Q , combined with a maximum force reso- 
lution of 2 ft, _1 kpc. 

• run B: baseline mass resolution unchanged, but the 
mass resolution of the galaxy models enhanced fur- 
ther by a factor of eight. 

• run C: same as above, but the gravitational soft- 
ening length for the particles in the galaxy models 
was decreased by a factor of two. 

• run D: baseline mass resolution enhanced further 
by a factor of eight to 1.50 x 10 7 h~ 1 MQ, mass 
resolution of the galaxy models equal to the base- 
line mass resolution, and a gravitational force res- 
olution equal to 1 /i _1 kpc everywhere in the high- 
resolution region. 

The runs A to D show progressively improving agree- 
ment as the resolution of the simulations is increased. 
This is especially evident for the smallest objects. The 
comparison of the runs C and D shows that the effects of 
two-body relaxation do not play an important role even 
in the lower resolution run C. The maximum resolution 
run exhibits a close agreement of the galaxy properties 
at z = with those assumed at the initial redshift. As in 
the case of Figure 1 , this figure shows a clear bend in the 
Kormendy relation but this time for the galaxy models 
that were not adiabatically compressed. 

In Figure 7, we show how the mass ratios evolve as a 
function of stellar mass. The color coding is the same as 
in Figure 1. The left panel is for the uncompressed initial 
models and the right one for the adiabatically contracted 
ones. As expected, the adiabatically contracted models 
have higher dark matter to stellar mass ratios. Neverthe- 
less, the conclusions remain the same. The obvious trend 
for the BCG to evolve in a dramatically different fashion 
than the rest of the population in the cluster is present. 
The conclusion regarding all other plots discussed above 
also hold for the models that were not adiabatically com- 
pressed. 

The mass of the cluster of galaxies considered here 
(10 15 M Q ) is much greater then the total mass in the 
entire computational volume (~4x 10 13 M Q ) considered 
in the cosmological simulations of the fundamental plane 
and scaling relations presented by Ohorbe et al. (2006). 
Their simulations include the effects of gas cooling and 
star formation. However, the mass resolution in our case 
is much better than the resolution of 2 x 64 3 particles 
these authors considered. It is not clear if their results 
arc numerically fully converged at z = as their con- 
vergence tests were performed at z = 1 due to the very 
high CPU requirements. Ohorbe et al. (2006) do not 
report a departure of the most massive object from the 
scaling relations but they do see a tilt in the fundamen- 
tal plane. In their model, the tilt is created during an 
early violent phase of rapid cooling and star formation. 



This tilted fundamental plane is then preserved in the 
dissipationlcss mergers that follow. 

We acknowledge that we can afford the very high nu- 
merical resolution and dynamical range we achieve in our 
cluster only at the expense of not including gas dynam- 
ical effects but, as a result, this allows us to address a 
very clean and well-defined problem and to isolate the 
purely dynamical effects from the complicated gas dy- 
namical processes. Also, in our approach we can ensure 
by construction that the galaxies at the initial high red- 
shift obey a prescribed mass-size relation, which further 
helps to make the problem we examine clean and well 
defined, and to arrive at robust conclusions. 

5. CONCLUSIONS AND DISCUSSION 

We have performed high dynamical range cosmological 
A-body simulation of the formation of a massive 10 15 M Q 
cluster of galaxies. We included both the stellar and dark 
matter particle components in order to investigate the 
(differential) evolution of these two components and the 
relationship between them. We also performed numerical 
convergence tests and checked the robustness of our ap- 
proach with respect to details of the initial galaxy mod- 
els by considering both regular and adiabatically com- 
pressed galaxy models. 

In our approach, the Kormendy relation (defined here 
as the relation between the half-light radius and the mass 
in the stellar component) evolves with redshift entirely 
due to stellar dynamical/ A-body processes. More specif- 
ically, the BCG clearly departs from the extrapolation of 
the initial relation to higher mass objects and the final 
relation appears to be bent at the massive end. The evo- 
lution of the BCG is accompanied by a systematic and 
strong increase in the dark matter-to-stellar mass ratio 
which is brought about by a significantly larger num- 
ber of mergers that the BCG experiences compared with 
other galaxies. In contrast, smaller galaxies show rather 
little evolution in this proxy of the mass-to-light ratio. 
We speculate that the above trend also implies that the 
slope of the R e — L relation for the BCGs is steeper than 
for the smaller early-type galaxies. This possibility would 
have to be explored in future simulations involving more 
clusters of varying mass. Interestingly, a gentle trend for 
this quantity to increase with the galaxy mass has been 
reported by Gallazzi et al. (2006) in early-type galaxies 
in the SDSS survey. The bend in the Kormendy relation 
has also been interpreted by Lauer et al. (2007) as being 
due to an increase in this mass ratio. We note that such 
variations may also have important implications for the 
lensing properties of the BCGs. 

Surprisingly, the Faber- Jackson relation only shows 
a tentative weak bending with redshift in our simula- 
tions that depends on the aperture size. However, iso- 
lated merger experiments performed mostly for parabolic 
merger orbits (Boylan-Kolchin et al. 2006) suggest that 
such an effect could be present. We speculate that the 
key factors responsible for this difference are: (1) the 
mergers in our model all take place in the central clus- 
ter potential and the BCG lies very close to it and (2) 
our initial conditions for the mergers are self-consistently 
taken into account as we evolve all galaxies in the cos- 
mological context starting from high redshift prior to the 
cluster assembly. We also note that the distribution of 
the initial orbital parameters of the merging galaxies is a 
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very sensitive function of eccentricity. More precisely, the 
distribution possesses a narrow spike at zero eccentricity 
(Benson 2005), and small departures from the idealized 
case of parabolic orbits are then likely in a realistic sce- 
nario, and (3) the precise definition of where the BCG 
envelope ends and the ICL component begins may some- 
what modify the measured strength of any bend in these 
relations. 
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Fig. 1. — Half-light radius vs. stellar mass relation. Filled red circles are for the initial models at z = 3 and the black ones are for z = 0. 
Intermediate redshift data is shown as open circles and distributed in redshift according to the color scale given in the upper left corner. 
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Fig. 2. — Mass in the dark matter (red) and stellar (violet) components for the BCG as a function of redshift. The values are normalized 
to the dark matter mass within the half-light radius at z = 3. 
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Fig. 3. — Cumulative number of mergers experienced by the three most massive galaxies in the cluster as a function of redshift. 
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Fig. 4. — Left panel: stellar mass as a function of the flux-weighted line-of-sight stellar velocity dispersion within the half-light radii. 
The meaning of symbols is the same as in Figure 1. Right panel: flux-weighted line-of-sight stellar velocity dispersion computed for three 
different apertures (half-light radius, 1/8 and 1/16 of that radius) as a function of the stellar mass at z = 0. 
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Fig. 5. — Fundamental "plane" for the galaxy model at z = 3 (red) and z = (violet). R e and cr e are the half-light radius and the 
velocity dispersion defined as in Figure 4. The intensity is defined as I e = m s tar/(7rijg). 
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Fig. 6. — Convergence study for the relation between half-light radius and stellar mass, R c — M*, for four different runs of increasing 
numerical resolution (from run A to Run D; see text for details). All filled circles correspond to the galaxies at z = 0. Open red circles are 
for the initial galaxy models at z = 3. All runs were carried out for uncontracted galaxy models. 
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Fig. 7. — Dark matter to stellar mass ratios vs. stellar mass. The meaning of the symbols is the same as in Figure 1. The left panel 
shows results for uncontracted galaxy models, while the panel on the right gives results for adiabatically compressed initial galaxy models. 



